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Based on first-principles density-functional theory calculations, we present a comparative study of the elec¬ 
tronic structures of ultranarrow zigzag graphene nanoribbons (ZGNRs) embedded in hexagonal boron nitride 
(BN) sheet and fully hydrogenated graphene (graphane) as a function of their width N (the number of zigzag 
C chains composing the ZGNRs). We find that ZGNRs/BN have the nonmagnetic ground state except at N 
= 5 and 6 with a weakly stabilized half-semimetallic state, whereas ZGNRs/graphane with N > 2 exhibit a 
strong antiferromagnetic coupling between ferromagnetically ordered edge states on each edge. It is revealed 
that the disparate magnetic properties of the two classes of ZGNRs are attributed to the contrasting interedge 
superexchange interactions arising from different interface structures: i.e., the asymmetric interface structure 
of ZGNRs/BN gives a relatively short-range and weak superexchange interaction between the two inequivalent 
edge states, while the symmetric interface structure of ZGNRs/graphane gives a long-range, strong interedge 
superexchange interaction. 

PACS numbers: 73.21.Hb, 73.22.Pr, 75.75.-c 


Graphene nanoribbons (GNRs) have been regarded as one 
of the most important classes of carbon-based nanomaterials 
due to their unique electronic and magnetic properties Fab¬ 
rication of GNRs with different widths and edges has been 
achieved by lithographic patterning^^ bottom-up fabrication^^ 
and chemical unzipping of carbon nanotubes It is known 
that the electronic and magnetic properties of GNRs vary with 
the ribbon width and the edge geometrythereby being uti¬ 
lized to design novel electronic and spintronic devices. For in¬ 
stance, the band gap of GNRs depends on the ribbon width^ 
and the zigzag-edged graphene nanoribbons (ZGNRs) have 
peculiar localized electronic states at both edges while the 
GNRs with armchair edges do not have such localized edge 
statesHere, the localized edge states of ZGNRs are ferro¬ 
magnetically ordered at each edge with an opposite spin orien¬ 
tation, forming an antiferromagnetic (AFM) spin order. Inter¬ 
estingly, it was predicted that such an AFM ordered ZGNR 
can have a half-metallic property if in-plane homogeneous 
electric field is applied across the edges of the ZGNR.- How¬ 
ever, the applied electric field is practically too high to re¬ 
alize half-metallic ZGNRsand therefore various alterna¬ 
tive approachesii"— have been proposed. Most of the alter¬ 
natives focus essentially on the same conceptual basis that the 
half-metallicity of ZGNRs can be enabled by the modifica¬ 
tion of edge states: e.g., the edge modification of ZGNRs 
with two different functional groups can produce the half- 
metallic property even in the absence of an external electric 
field^i^iii Such an asymmetric edge modification can also be 
achieved when ZGNRs are embedded in a hexagonal boron 
nitride (/z-BN) sheet These embedded ZGNRs (hereafter 
designated as ZGNRs/BN) have a C—B interface at one edge 
and a C—N interface at the opposite edge [see Fig. 1(a)]. A 
number of density functional theory (DFT) calculations re¬ 


ported the presence of half-metallicity or half-semimetallicity 
in ZGNRs/BN On the experimental side, the fabrication 
of graphene-BN hybrid structures was recently reported, with 
graphene strips as narrow as tens of nanometers In a dif¬ 
ferent way, the ultranarrow ZGNRs can be fabricated2£ by re¬ 
moving hydrogen atoms from a fully hydrogenated graphene 
(viz. graphane)}^! see Fig. 1(b). Such embedded ZGNRs 
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FIG. I: (Color online) Top view and side view of the optimized struc¬ 
tures of (a) the ZGNR/BN and (b) the ZGNR/graphane with N = 5. 
The V and y axes are taken to be parallel and perpendicular to the 
edges of ZGNR, respectively. L and R represent the left and right 
edges, respectively. 



Typeset by REVT^X 


2 


> 

U) 

Q> 

C 

lU 


> 

3 - 

>* 

a> 

£ 

lU 




FIG. 2: (Color online) Calculated band structures of (a) ZGNRs/BN and (b) ZGNRs/graphane as a function of N. The results for the h-BN 
sheet and graphane are also given. The charge character of edge states is shown with an isosurface of 0.01 electrons/A^. In the schematic 
diagram, Cb (Cn) indicates the edge C atom bonding to a B (N) edge atom. The energy zero represents the Fermi level. The direction of F — 
X line is parallel to the edges. For distinction, C atoms composing the ZGNR with = 1 are drawn with circles in different brightness. 


(hereafter designated as ZGNRs/graphane) were predicted to 
exhibit the insulating AFM ground statcr^r^i similar to iso¬ 
lated ZGNRsii Note that ZGNRs/graphane have the symmet¬ 
ric interface structure with an identical C—CH interface on 
both edges, differing from the asymmetric interface structure 
of ZGNRs/BN. Therefore, it is very interesting to explore the 
roles of interface structure in determining the drastically dif¬ 
ferent electronic and magnetic properties of ZGNRs/BN and 
ZGNRs/graphane. 

In this paper, we perform first-principles DFT calcula¬ 
tions to investigate the electronic structures of ultranarrow 
ZGNRs/BN and ZGNRs/graphane as a function of N from 
1 to 8. We find that, in contrast to ZGNRs/graphane where 
the nonmagnetic (NM) ground state aX N = 1 is converted to 
the AFM ground state for N > 2, ZGNRs/BN have the NM 
ground state except at A = 5 and 6 with the half-semimetallic 
state. Such different behaviors of the ZGNRs embedded in 
h-BN sheet and graphane can be traced to the contrasting 
features of edge states due to their different interface struc¬ 
tures. The asymmetric interface structure of ZGNRs/BN pro¬ 
duces the inequivalent edge states originating from the B—C 
and N—C interfaces, giving rise to a relatively short-range 
and weak superexchange interaction between the two edge 
states. On the other hand, the symmetric interface structure 
of ZGNRs/graphane produces the identical edge states with 
partially fiat bands, leading not only to a magnetic instability 
due to the enhanced density of states (DOS) at the Fermi level 
(Ef) but also to a long-range, strong interedge superexchange 
interaction. 

The present first-principles DFT calculations were per¬ 


formed using the Fritz-Haber-Institute ab-initio molecular 
simulations (FHI-aims) code^^ for an accurate, all-electron 
description based on numeric atom-centered orbitals, with 
“tight” computational settings and accurate tier-2 basis sets. 
For the exchange-correlation energy, we employed the gen¬ 
eralized gradient approximation functional of Perdew-Burke- 
Ernzerhofi2& The embedded ZGNRs were simulated using a 
periodic supercell with a constant in-plane unit cell length of 
~46 A (changing the width of BN or graphane from ^26 
to ^41 A with respect to N) and a vacuum spacing of ^30 
A between the periodic sheets. For the Brillouin zone inte¬ 
gration, we used 128x1x1 k-points in the surface Brillouin 
zone. All of the atoms were allowed to relax along the cal¬ 
culated forces until all the residual force components are less 
than 0.02 eV/A. The optimized lattice constant and band gap 
of h-BN sheet (graphane) are found to be 2.51 (2.54) A and 
4.66 (3.62) eV, respectively, in good agreement with previous 
DFT calculations 

We begin to determine the atomic and electronic structures 
of ZGNRs/BN and ZGNRs/graphane using spin-unpolarized 
calculations. The calculated band structures of ZGNRs/BN 
and ZGNRS/graphane are displayed as a function of N in 
Fig. 2(a) and 2(b), respectively, together with those of h- 
BN sheet and graphane. For the ZGNR/graphane with N = 
1, we obtain a bond-alternated structure with two different 
C—C bond lengths, dc-c = 1-42 and 1.45 A [see the inset 
of Fig. 2(b)], indicating a Peierls distortion of Ad = ±0.015 
A. This Peierls distortion accompanies a band-gap opening of 
0.22 eV between the 7Z and tt* bands. On the other hand, for 
the ZGNR/BN with A = 1, such a bond alternation does not 
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TABLE I: Calculated energy difference (in meV/unit cell) between the NM and half-semimetallic (HS) or AFM configurations for ZGNRs/BN 
and ZGNRs/graphane as a function of N. The band gap in each system is also given in the unit of eV. 
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ZGNRs/BN 

ground state 

NM 

NM 

NM 

NM 

HS 

HS 

NM 

NM 


AFnm-hs 

- 

- 

- 

- 

1.7 

1.5 

- 

- 


Eg 

1.711 

0.701 

0.176 

0.020 

0.023 

0.002 

0.005 

0.010 

ZGNRs/graphane 

ground state 

NM 

AFM 

AFM 

AFM 

AFM 

AFM 

AFM 

AFM 


AFnm-afm 

- 

8.0 

32.6 

50.3 

59.5 

64.5 

69.1 

72.7 


Eg 

0.217 

0.410 

0.532 

0.533 

0.509 

0.479 

0.450 

0.422 


occur with an equal C—C bond length of dc-c = 1-43 A, and 
the charge character of the n (tt*) state at the X point repre¬ 
sents the hybridization between C and B (N) atoms [see the 
inset of Fig. 2(a)], giving rise to a large band-gap opening of 
1.71 eV. As N increases, the C—C bond lengths in ZGNRs/BN 
and ZGNRs/graphane are close to each other as 1.44 and 1.45 
A, respectively. Figure 2(a) shows that, as N of ZGNRs/BN 
increases, the band gap (Eg) decreases, almost being closed 
from = 4. However, we note that the breaking of the sublat¬ 
tice symmetry in ZGNRs/BN, due to their asymmetric inter¬ 
face structure, avoids the crossing of n and tt* bands at Ef 
[see the inset of A^ = 8 in Fig. 2(a)]. On the other hand, 
for ZGNRs/graphane with identical edge interfaces, the % and 
TT* bands cross the Fermi level with increasing N, formimg a 
partially two-fold degenerate flat band at a sufficiently wider 
width [see Fig. 2(b)]. 

It is noteworthy that the charge characters of 7Z and TT* states 
in ZGNRs/BN represent asymmetric edge states localized at 
the C—B and C—N interfaces, respectively [see the inset for 
A^ = 8 in Fig. 2(a)]. Here, the schematic diagram of fron¬ 
tier orbital interactions shows that the highest occupied 7Z and 
lowest unoccupied TT* states are characterized as the C—B 
bonding and C—N antibonding orbitals, respectively. On the 
other hand, the ZGNR/graphane with N = S shows symmet¬ 
ric TT and TT* edge states localized at both sides [see the in¬ 
set for A^ = 8 in Fig. 2(b)]. These different features of edge 
states between ZGNRs/BN and ZGNRs/graphane may influ¬ 
ence the range and strength of the interaction between two 
edges. In order to compare the effects of the interedge in¬ 
teraction on the half-semimetallicity of ZGNRs/BN and the 
AFM order of ZGNRs/graphane, we perform spin-polarized 
calculations for the two systems as a function of N. It is 
known that the electric fleld created by different electrostatic 
potentials at the C—B and C—N interfaces is associated with 
half-semimetallicity in ZGNRs/BNr^ while the flat-band-like 
character in the edge states of ZGNRs/graphane induces a 
magnetic instability due to the enhanced DOS at Ep^ The 
calculated stabilization energies of half-semimetallicity and 
AFM order in ZGNRs/BN and ZGNRS/graphane relative to 
the corresponding NM conflguration are given as a function 
of N in Table I. In ZGNRs/BN, we obtain the NM ground 
state for A^ < 4 and A^ > 7, while the half-semimetallic ground 
state at A^ = 5 and 6. This trend showing that wide-(or ex¬ 
tremely narrow) and intermediate-width ZGNRs are stabilized 
as the NM and half-semimetallic conflgurations, respectively. 


is consistent with previous DFT studyJ^ Note that the half- 
semimetallic conflguration at A^ = 5 and 6 is only a few meV 
lower in energy than the corresponding NM conflguration 
(see Table I), indicating that the interedge interaction produc¬ 
ing half-semimetallicity in ZGNRs/BN is very weak. On the 
other hand, in ZGNRs/graphane, we obtain the AFM ground 
state for N >2, where the total-energy difference AE’nm-afm 
between the NM and AFM conflgurations monotonicahy in¬ 
creases as N increases, reaching ^73 meV at A^ = 8 (see Table 
I). These results obviously indicate that the interedge interac¬ 
tion in ZGNRs/graphane is long-range and strong compared to 
that in ZGNRs/BN. It is remarkable that the geometric sym- 


(a) 



FIG. 3: (Color online) Calculated AFM band structures of (a) 
ZGNR/BN and (b) ZGNR/graphane for A = 5. The energy zero 
represents the Fermi level. The spin densities of ZGNR/BN and 
ZGNR/graphane are also given. The inset in (a) magnifies the band 
gap of the spin-down bands. The spin densities are drawn with an 
isosurface of 0.02 (—0.02) electrons/ A^ . 
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metry of two edges in ZGNRs embedded in either /z-BN sheet 
or graphane plays an important role in determining the range 
and strength of interedge interaction. 

Figure 3(a) and 3(b) show the comparison of the 
spin-polarized band structures of the ZGNR/BN and 
ZGNR/graphane with N = 5. In the former system, the spin- 
up and spin-down bands open a gap of 0.14 and 0.02 eV 
[see Fig. 3(a)], respectively. These values of gap opening 
are much smaller compared to the ZGNR/graphane system 
where the spin-up and spin-down bands open an identical 
band gap of 0.51 eV [see Fig. 3(b)]. Here, the much smaller 
band gap in the half-semimetallic ZGNR/BN compared to the 
AFM ZGNR/graphane gives rise to a much smaller value of 
A£’nm-hs = 1.7 meV than AEnm-apm = 59.5 meV (see Table 
I). It is interesting to notice that there is a subtle difference of 
the spin characters between the half-semimetallic ZGNR/BN 
and the AFM ZGNR/graphane. As shown in Fig. 3(a) and 
3(b), the spin density of the former system is relatively well- 
localized around the two edges, whereas that of the latter sys¬ 
tem shows some extension up to the middle of the ribbon. 
This reflects relatively short-range (weak) versus long-range 
(strong) interedge spin-spin interactions in ZGNRs/BN and 
ZGNRs/graphane. 

To understand the microscopic mechanism for the 
half-semimetallicity and AFM order in ZGNRs/BN and 
ZGNRs/graphane, we plot, in Fig. 4(a) and 4(b), the spin- 
polarized local DOS projected onto the two edge C atoms [in 
the left (L) or right (R) edge site in Fig. 1] together with their 
spin characters. For the ZGNR/BN with A = 5, it is seen that 
the occupied (unoccupied) spin-up and spin-down edge states 
are localized at the L (R) edge [see Fig. 4(a)]. However, for 
the ZGNR/graphane with A = 5, the occupied (unoccupied) 
spin-up and spin-down edge states are localized at the L (R) 
and R (L) edges, respectively [see Fig. 4(b)]. Since electronic 
states with the same spin direction can hybridize with each 
other, the hybridization occurs between the occupied and un¬ 
occupied spin-up or spin-down states localized at the L and 
R (or R and L) edges. This kind of exchange interaction be¬ 
tween the occupied and unoccupied states is characterized as 
a superexchange mechanism.i^L^ Such an interedge superex¬ 
change interaction leads to a relatively long-range, strong 
interedge spin-spin interaction in ZGNRs/graphane, thereby 
giving rise to a large energy gain in AEnm-apm as well as a 
large gap opening, as shown in Table I. 

In summary, using first-principles DFT calculations, we 
have performed a comparative study of the electronic struc¬ 
tures of ultranarrow ZGNRs/BN and ZGNRs/graphane. Such 
embedded ZGNRs in /z-BN sheet and graphane are found 
to exhibit drastically different electronic characteristics. Un¬ 
like ZGNRs/graphane, whose NM configuration exhibits par¬ 
tially flat bands at Ef as A increases, ZGNRs/BN do not 
have such a flat-band-like character. Consequently, the for¬ 
mer ZGNRs show a magnetic instability due to the enhanced 


DOS at Ef, whereas the latter ZGNRs preserve the NM 
ground state except at A = 5 and 6 with a half-semimetallic 
state. We revealed that the disparate magnetic properties of 
the two classes of ZGNRs can be traced to the different fea¬ 
tures of their interface structures: i.e., unlike the symmetric 
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FIG. 4: (Color online) The spin-polarized local DOS projected onto 
the two edge C atoms [in the left (L) or right (R) edge site in Fig. 1] 
of (a) ZGNR/BN and (b) ZGNR/graphane with N = 5. The energy 
zero represents Fermi level. The charge characters of the spin-up and 
spin-down states for the occupied and the unoccupied band are taken 
at the X point with an isosurface of 0.04 (—0.04) electrons/ 


interface structure of ZGNRs/graphane, the asymmetric inter¬ 
face structure of ZGNRs/BN produces the inequivalent edge 
states on both sides of the nanoribbons, giving rise to a rela¬ 
tively short-range, weak interedge superexchange interaction. 
The resulting different electronic and magnetic properties of 
ZGNRs/BN and ZGNRs/graphane may be utilized for the ap¬ 
plication of nano-scale electronic devices such as conducting 
wires or field-effect transistors. 
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